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Abstract. In dissipative dynamical systems phase space volumes contract, on average. 
Therefore, the invariant measure on the attractor is singular with respect to the Lebesgue 
measure. As noted by Ruelle, a generic perturbation pushes the state out of the attractor, 
hence the statistical features of the perturbation and, in particular, of the relaxation, cannot 
be understood solely in terms of the unperturbed dynamics on the attractor. This remark 
seems to seriously limit the applicability of the standard fluctuation dissipation procedure in 
the statistical mechanics of nonequilibrium (dissipative) systems. In this paper we show that 
the singular character of the steady state does not constitute a serious limitation in the case of 
systems with many degrees of freedom. The reason is that one typically deals with projected 
dynamics, and these are associated with regular probability distributions in the corresponding 
lower dimensional spaces. 
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1. Introduction 

Since its early developments, due mainly to the works of L. Onsager and R. Kubo [HI El El SJ, 
the fluctuation-dissipation theorem (FDT) represents a cornerstone in the construction of 
a theory of nonequilibrium phenomena f5|. This celebrated result was developed in the 
context of Hamiltonian dynamical systems, slightly perturbed out of their thermodynamic 
equilibrium, and it was later extended to stochastic systems obeying a Langevin Equation 
lis The importance of the FDT rests on the fact that it sheds light on the crucial relation 
between the response Rv{t) of a system to an external perturbation and a time correlation 
function computed at equilibrium. In other words, having perturbed a given Hamiltonian 
Hq with an external field h^, to obtain the perturbed Hamiltonian Hq — heV, where V is an 
observable conjugated with he, the FDT allows us to compute nonequilibrium quantities, such 
as the transport coefficients [I8l|9l|, solely in terms of the unperturbed equilibrium state. On the 
other hand, a generic dynamical system is not Hamiltonian: for phenomenological practical 
pruposes, one typically deals with dissipative dynamics, as in the important case of viscous 
hydrodynamics 

The invariant measure of a chaotic dissipative system, /i say, is singular with respect 
to the Lebesgue measure and is usually supported on a fractal attractor. This is not just a 
mathematical curiosity, it is a potential source of difficulties for the applicability of the FDT 
in dissipative systems. Indeed, the standard FDT ensures that the statistical features of a 
perturbation are related to the statistical properties of the unperturbed system, but that cannot 
be the case in general, in dissipative systems. The reason is that, given an initial state x{0) on 
the attractor and a generic perturbation 5x{0), the perturbed initial state x^^^ (0) = x(0) +5x{0) 
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and its time evolution may lie outside the support of /i, hence their statistical properties 

cannot be expressed by fi, which attributes vanishing probability to such states. In the 

cases considered by Ruelle MOI . the perturbation (5x(0) and its time evolution 5x_{t) can 

be decomposed as the sum of two parts, 5x^{t) and respectively perpendicular and 

parallel to the "fibres" of the attractor, 

5x{t) = Sx^{t) + 

which makes it natural to expect that the statistical features of be related to the 

dynamics on the attractor, while it is easy to construct examples in which Sxj_{t) is not 
described by the unperturbed dynamics. 

From the mathematical point of view, this fact is rather transparent. On the other hand, it 
should not be a concern in statistical mechanics, except in pathological cases. Indeed, a series 
of numerical investigations of chaotic dissipative systems shows that the standard FDT holds 
under rather general conditions, mainly if the invariant measure is absolutely continuous with 
respect to Lebesgue, cf. Ref.[6l| for a review. Moreover, although dissipative systems have 
singular invariant measures, any small amount of noise produces smooth invariant measures, 
which allow generalized FDTs to be expressed solely in terms of the unperturbed states, 
analogously to the standard equilibrium case. Apart from technical aspects, the intuitive 
reason for which the FDT in systems with noise can be expressed only in terms of the invariant 
measure, is that 3:^(0) remains within the support of this measure. 

In this paper, we want to take advantage of the fact that a similar situation is realized 
without any noise, if one works in the projected space of the physically relevant observables. 
Indeed, marginals of singular phase space measures, on spaces of sufficiently lower dimension 
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than the phase space, are usually regular lfTTl[T2]| . 

Our paper is organized as follows: Section |2] is devoted to a short presentation of some 
general results on FDT for chaotic dissipative systems. In Sec. [3] we discuss the numerical 
results for two dissipative chaotic maps, showing that the singular character of their invariant 
measures does not prevent the response of standard observables to be expressed only in terms 
of the invariant measure, as in the standard case. Conclusions are drawn in Sec. IH 

2. Some results on FDT in chaotic dissipative systems 

Let us concisely recall Ruelle's approach to linear response in deterministic dissipative 
dynamical systems ifTOll . Let {M, 5*, /i) be a dynamical system, with M its compact phase 
space, S** : — 7- a one parameter group of diffeomorphisms and /i the invariant natural 
measure. Following Ruelle [lOJ, who considers axiom A systems, one may show that the 
effect of a perturbation = 5F\\^{t)+5Fx_{t) on the response of a generic (smooth enough) 
observable A attains the form: 



where the subscript | refers to the dynamics on the unstable tangent bundle (along the 
attractor), while _l refers to the transversal directions, cf. left panel of Fig.[IJ Ruelle's central 

{A) 

remark is that i^y may be expressed in terms of a correlation function evaluated with respect 
to the unperturbed dynamics, while R^f^ depends on the dynamics along the stable manifold, 
hence it may not be determined by fi, and should be quite difficult to compute numerically 




(1) 



m. 
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Figure 1. Left panel: In Ruelle's approach, the perturbation is expressed as the sum of 
one component parallel to the unstable manifold and one parallel to the stable manifold. 
Right panel: In the present work, the reference frame is rotated so that the direction of the 
perturbation coincides with one of the basis vectors. 

To illustrate these facts, the Authors of Ref. ffT3l study a 2-dimensional model, which 
consists of a chaotic rotator on a plane and, for such a system, succeed to numerically 
estimate the R^f^ term in eq.©. Nevertheless, in the next Section, we argue that R^f^ may 
spoil the generalized FDT only if the perturbation is carefully oriented with respect to the 
stable and unstable manifolds. This is only possible in peculiar situations, such as those of 
Ref. [fT3l . in which the invariant measure is the product of a radial and and angular component 
and, furthermore, the perturbation lies on the radial direction, leaving the angular dynamics 
unaffected. 

A different approach to the FDT has been proposed in lfT4l . which concerns deterministic 
dynamics perturbed by stochastic contributions. Here, the invariant measure /i can be assumed 
to have density p: d/i(x) = p{x)dx. Then, if the initial conditions are modified by an 
impulsive perturbation a^Q — Xq + 5xq, the invariant density p(xq) is replaced by a perturbed 
initial density Po{xq]5xq) = p{xq — 5xq), where the subscript denotes the initial state, 



Fluctuation-dissipation relation for chaotic non-Hamiltonian systems 6 
right after the perturbation. This state is not stationary and evolves in time, producing time 

dependent densities pt{xi^', 5"^), which are assumed to eventually relax back to p{xq). Given 

the transition probability W{xj^,Q — )■ t) determined by the dynamics, the response of 

coordinate Xi is expressed by: 

5xi{t) ^ j j [pfeo ~ ^^o) ~ P(^o)] ^feo' — )■ X, t)dxQdx (2) 
and one may introduce the response function Rij as [[T4| : 

>■ n \ 

(3) 



t=0 



dxj 

which is a correlation function computed with respect to the unperturbed state. It is worth 
to note that it makes no difference in the derivation of eq.© whether the steady state is an 
equilibrium state or not; it suffices that p be differentiable. 

Let us consider again Eq.© and, for sake of simplicity, assume that all components of 
Sx{0) vanish, except the i-th component. Then, the response of Xi may also be written as: 

XiBi{xi,Sx^,t)dxi (4) 
where Bi{xi, 5xq, t), defined by the term within curly brackets, may also be written as: 

Bi{xi, 5xj^, t) = pt{xi; 5xq) - p{xi) (5) 
where p{xi) and pt are the marginal probability distributions defined by: 

As projected singular measures are expected to be smooth, especially if the dimension 
of the projected space is sensibly smaller than that of the original space, one may adopt the 
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same procedure also for dissipative deterministic dynamical systems. Indeed, the response 

function Bi{xi, Sxq, t) in Eq.© is also expected to be smooth, and to make the response of Xi 

computable from the invariant measure only. In the next section we investigate this possibility. 

3. Coarse graining analysis 

In terms of phase space probability measures, the response formula Eq.© reads: 

5xi{t) = / d/ii(x; 5xo) - / d/i(x) (6) 



where d/it(x; 5xq) is the time evolving perturbed measure whose initial state is given by 

Because dissipative dynamical systems do not have an invariant probability density, it is 
convenient to introduce a coarse graining in phase space, to approximate the singular invariant 
measure /i by means of piecewise constant distributions. 

Let us consider a ci-dimensional phase space M, with an e-partition made of a finite set 
of (i-dimensional hypercubes Afc(e) of side e and centers Xf,. Introduce the e-coarse graining 
of /i and of Ht defined by the probabilities Pki^) and Pt,k{^] 5xq) of the hypercubes Afc(e): 

Pfc(e) = / d/i(x) , Pt,fc(e;5^o) = / d/it(x; feg) . (7) 
This leads to the coarse grained invariant density e): 

Pfc(e)/e'^ ifxeAfc(e) 
Pfe = 2^ Pfcfe e) , with pfc(x; e) = <j (8) 

^ I else 

Let Zi be the number of bins of of form xf^ — e/2, xf'' + e/2 j , g G {1,2,..., Zj}, in the z-th 
direction. Then, the marginalization of the coarse grained distribution yields the following set 



Fluctuation-dissipation relation for chaotic non-Hamiltonian systems 
of Zi probabilities: 



p-''^(e) = y^^^ I J p{x; e) Y[dxj > dxi = Prob (^xi G 



4'' -'^,4'^+'^)) (9) 



each of which is the invariant probability that the coordinate Xi lie in one of the Zi bins. 
In an analogous way, one may define the marginal of the evolving coarse grained perturbed 
probability p[f (e; 6xq). In both cases, dividing by e, one obtains the coarse grained marginal 
probability densities pf\e) and pf'^ {e; 6xq), as well as the e-coarse grained version of the 
response function Bi{xi, 6xj),t): 



Bl'^\x,,6x„t,e) = ^ Pl'::{e,6x,)-p^>{e) = pfj{e,5x,)-pf\e) (10) 



1 



In the following, we will show that the r.h.s. of Eq.ClOl) tends to a regular function of Xj in the 
— oo, e — 0, limit. Then, in the limit of small perturbations feg, b\'^^ (xj, feg, t, e) may be 
expanded as a Taylor series, to yield an expression similar to standard response theory, in the 
sense that it depends solely on the unperturbed state. The difference, here, is that the invariant 
measure is singular and represents a nonequilibrium steady state. 

To illustrate this fact, we run a set of trajectories with uniformly distributed initial 
conditions in the phase spaces of two simple, but substantially different, 2-dimensional maps: 
a dissipative baker map, and the Henon map. 
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3.1. The dissipative baker map 



Let M = [0,1] X [0, 1] be the phase space, and consider the evolution equation 

for < Xn < I', 



I \ 

y 2/n.+ l j 



M 



( \ 



whose Jacobian determinant is given by 

J A = r/l 

Juix) 



\ ryn J 

r + IVn 



for < a: < /; 



.(11) 



for I < Xn < \. 



(12) 



= Ijr = J^S for / < a; < 1. 
and shows that the M is dissiaptive for / < 1/2. The map M is hyperbolic, since stable 

and unstable manifolds which intersect each other orthogonally are defined at all points 

X G A^, except in the irrelevant vertical segment at x = /. The directions of these manifolds 

coincide, respectively, with the vertical and horizontal directions. It can also be shown that 

this dynamical system is endowed with an invariant measure \i which is smooth along the 

unstable manifold and singular along the stable one, cf. FigsjU In particular, /i factorizes as 

d/i(x) = dx X dA(y), similarly to the case of lfT3l . 

In order to verify whether the functions corresponding to the above introduced 

(5xq, t, e) become regular functions in the fine graining limit, let us consider first 

an impulsive perturbation, directed purely along the stable manifold, i.e. bx^ = {0,5yo). 

Ruelle's work on singular measures is clearly relevant, in this case, because the support of the 

marginal perturbed probability measure, obtained projecting out the y-direction has simply 

drifted preserving its singular character, while the state may have fallen outside the support of 

the unperturbed invariant measure, cf. left panel of Fig (31 
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Figure 2. Invariant probability distribution of the map defined by Eq. (fTTl l. 

Consider now an initial impulsive perturbation with one component, no matter how 
small, along the unstable manifold, 5xq = (5xo, 5?/o) and rotate the vectors of the basis of 
the 2-dimensional plane, so that the coordinate x lies along the direction of the perturbation, 
as shown in the right panel of FiglU We find that Bx\x, 6xj^,t, e) is regular as a function 
of X. Indeed, the projections of /i and of its perturbations onto the direction of 6xq have a 
density along all directions except the vertical one, cf. right panel of FigjSl Hence, a small 
perturbation does not take the state outside the corresponding projected support. 

As already noted in lfT3l . this Baker map shows that the response to very carefully 
selected perturbations, cannot be computed in general from solely the invariant measure. 
However, similarly to the case of lfT3l . the factorization of /i makes the present case rather 
peculiar. Indeed, for the overwhelming majority of dynamical systems, it looks impossible 
to select directions such that the projected measures preserve the same degree of singularity 
as the full measures. This is a consequence of the fact that stable and unstable manifolds 
have different orientations in different parts of the phase space, provided they exist. Clearly, 
the higher the dimensionality of the phase space and the larger the number of projected out 
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Figure 3. Left panel: singular distribution of the unperturbed system, obtained projecting the 
invariant measure onto the vertical direction y, which is the direction of the stable manifold. 
Right panel: projected invariant density along the direction of the perturbation, which forms 
an angle a = 7r/8 radiants with the y-direction. z = 1 — (a;cos(a) — 2/sin(a) + sin(a)) is 
the coordinate along this direction. Although hardly differentiable, this projected distribution 
has a density. 

dimensions, the more difficult it is to preserve singular characters. 
3.2. The Henon map 



Consider for instance the Henon map defined by: 

/ \ / 

M 

\ Vn+l j 



( \ 



\ 



i/„ + 1 - ax\ 



\ 
I 



(13) 



one the phase space M. = [— |, |] x [|, |], where a = 1.4 and b = 0.3 imply a chaotic 
dissipative dynamics, with a fractal invariant measure fi, which is not the product of the 
marginal measures obtained by projecting onto the horizontal and the vertical directions. 
These marginals are indeed regular and would yield a regular product. As stable and 
unstable manifolds wind around, changing orientation, in a very complicated fashion, it seems 



impossible, here, to disentangle the contributions of one phase space direction from the other. 
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Figure 4. Left panel: plot of the function Bx [x, 5x^,1, e), introduced by Eq.(|5]l, for 
the Henon map, after t = 2 iterations, with an initial perturbation along the x direction, 
^Sio = (2-5 • 10^^,0), and different numbers of bins Zx on the a;-axis, as well as different 
number N of trajectories, so to keep the statistics in the e — > limit: = 3 • 10®, Zx = 1320 
(green curve), iV = 5 • 10®, Zx = 1650 (purple curve), 8 • 10^ Zx = 1980 (blue curve). 
Right panel: plot of the function B^,'''' {y, 6x^,1, e) after t = 2 iterations, for the same initial 
perturbation considered in the left panel, and different values of bins Zy on the ?;-axis and 
different values of iV: iV = 3 • 10^ Zy = 336 (green curve), iV = 5 • 10^ Zy = 420 (purple 
curve), iV = 8 ■ 10®, Zy = 505 (blue curve). The curves largely overlap, but the figure does 
not clarify whether si"^' , By"^^ get smoother as Zx and Zy increase with . 

Then, because no direction appears to be priviledged in phase space, an initial 
perturbation along one of the axis should not lead to any singular perturbed projected measure, 
or irregular response function, see e.g. Refs. [fTn[T2ll . Unfortunately, this is not obvious from 
the histograms constructed with growing numbers of bins, as they seem to be quite irregular 
and to develop singularities in some parts of the phase space, cf. Figs. |4] and [5l However, this 
does not necessarily prevent the projected measures from having a density. 
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Figure 5. Plot of the projected invariant probability density (e), for the Henon map, with 
N = b- 10^ = 1650 (purple curve, with green errorbars) and TV = 8 • 10^ = 1980 
(blue curve). The figure does not clarify whether the invariant distribution is siungular or not. 

Therefore, to clarify whether the projected probability density pi?^ (e) exists or not in the 
e — limit, we have examined the behavior of the Shannon Entropy, defined as 

S^{e) = -eY,p'f\e)\og{p^\e)) (14) 

q=l 

with e the size of the bin along the direction of the perturbation. Note that this entropy is often 
defined differently; our definition is meant to introduce a quantity whose e — )■ limit is finite 
if a density exists, while it diverges if the measure is singular. We approximated Sx{^) by 
running different sets of trajectories, with different sizes of the coarse graining of the x-axis. 
Our simulations with = 2 ■ 10^, show that Sx has substantially converged to its asymptotic 
— 7- oo limit, cf. Fig. [6l Moreover, for fixed N , Sx decreases as the number Zx of bins 
grows, and appears to tend to a constant as l/Z^. — )■ 0, cf. Fig.|7l 

Figures [8] and [9] further prove that Sa is always a finite quantity in the Henon case while, 
in the baker case, it diverges logarithmically only when Sa tends to n/2, which is the only 
angle for which the projected invariant measure is singular. Therefore, the response can be 
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0.845 




Figure 6. Numerical simulations of the Shannon Entropy S'x(e) for the Henon map, where 
denotes the number of bins considered on the x-axis and N denotes the number of trajectories. 
The curves collapse onto each other, approximating the asymptotic value of Sx associated with 
the N, Zx ^ oo limits. 




Figure 7. Left panel: Fit of the numerical data for Sx corresponding to = 7 • 10^, with 
f{x) = flo + "l^, ao = 0.7998 and ai ~ 5.3847. Right panel: Same curve as in the left panel, 
plotted vs. the variable e — l/Zx, to extrapolate the asymptotic value Sx — 0.8. 
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Figure 8. Numerical data with error bars, for the Shannon Entropy Sa of the Henon map ( fT3] ). 
as a function of the angle a E [0, 7r/2] with respect to the horizontal direction, at increasing 
orders of magnification. The top right panel is a magnification of the framed part of the top left 
panel. The bottom left panel is a magnification of the framed part of the top right panel. The 
bottom right panel is a magnification of the framed part of the bottom left panel. The quantity 
Sa is quite structured, especially for a close to 7r/2, but it is a finite quantity. 

obtained from the invariant measure at all perturbation angles in the case of the Henon map, 
and at all but a single angle for the Baker map. This confirms the applicability of a generalized 
FDT, which yields the response function in terms of the unperturbed state only, even if 
supported on a fractal set, except in very special situations, such as a negligible set in cases in 
which the invariant measure is the product of regular and singular mesaures. In particular, for 
the baker map, the response to a perturbation may be expressed just in terms of the smooth 
projected invariant measure if one does not perturb uniquely the vertical coordinates. For the 



Henon map, all directions lead to the existence of a projected invariant measure, although 
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Figure 9. Left panel: Numerical data, supplemented by error bars, for the Shannon Entropy 
Sa of the Baker map ( fTTT i. as a function of the angle a e [0, 7r/2] with respect to the x- 
axis. Right panel: Fit of for the Baker map, with the curve /(a) = alog(^ — a) + b, 
a = -0.101524 ± 0.002768, and b = 0.0411212 ± 0.01401. In this case, is a much 
simpler function of a than in the Henon case, and it diverges in the a ^ Tr/2 limit, because 
the projected invariant measure is singular at a = 7r/2. 



very finely structured. 



4. Conclusions 

In this paper we liave reviewed the methods proposed by Ref. [[TOl and by Refs. [[61 [T4l . 
concerning the derivation of response formulae for systems in nonequilibrium steady states. 
In particular, we have shown that the idea of [[T4l . which is based on the existence of a smooth 
invariant probability density, may be applied quite generally, with suitable adjustments, to 
dissipative deterministic dynamics. This requires that projected distributions be considered, 
rather than the full phase space distributions, becuase projected distributions are usually 
regular. Only very special combinations of dynamics and perturbations seem to prevent 
this approach, although in low dimensional dynamics such as ours, the projected distribution 
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functions appear to be quite complex and not smooth. 

The presence of noise, in any physically relevant dynamical system, does contribute 
to smooth out the invariant density, but even in the absence of noise, the fact that statistical 
mechanics is typically interested in projected dynamics allows an approach to FDT which only 
requires the properties of the unperturbed states, as in standard response theory. Clearly, this 
is better and better justified as the dimensionality of the phase space grows. In particular, it is 
appropriate for macroscopic systems in nonequilibrium steady states, because the dynamics 
of interest take place in a space whose dimensionality is enormously smaller than that of 
the phase space. Then, as projecting out more and more produces smoother and smoother 
distributions, one finds that the approach of Ref.[[T4| can be used to obtain the linear response 
function about nonequilirium steady states, from the unperturbed measure only. 

Our results support the idea that the projection procedure makes unnecessary the explicit 
calculation of the term discovered by Ruelle, which was supposed to forbid the standard 
approach. This does not mean that Ruelle's term is necessarily negligible [13]. However, 
except in very peculiar situations, such as our maker map which has carefully oriented 
manifolds, and for carefully chosen perturbations, that term does not need to be explicitly 
computed and the calculation of response may be carried out referring only to the unperturbed 
dynamics, as in the standard cases. 
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